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ABSTRACT 

We consider a simplest extension of the Standard Model (SM) through the incorporation of a 
real scalar singlet and an additional discrete Z 2 symmetry. The model admits the neutral scalar 
singlet to be stable and thus, a viable component of dark matter. We explore the parameter space 
of the model keeping in view the constraints arise from different dark matter direct detection 
experiments through WIMP-nucleon scattering. First of all, we have utilised the data obtained 
from CDMS, XENON- 10 and XENON- 100 collaborations. We further constraint the parameter 
space from the DAMA collaboration results (both with and without channelling) and CoGeNT 
collaboration results. Throughout our analysis, the constraint arises due to the observed relic 
density of dark matter reported by WMAP experiment, is also incorporated. Utilising all those 
constraints, on the model parameter space, we calculate the event rates and the annual variation 
of event rates in the context of a Liquid Argon Detector experiment. 
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1 Introduction 



Several cosmological observations like rotation curves of spiral galaxies, the gravitational micro- 
lensing, observations on Virgo and Coma clusters [1, 2], bullet clusters [3], etc. provide indications 
of existence of huge amount of non-luminous matter or dark matter(DM) in the universe. The 
Wilkinson Microwave Anisotropy Probe (WMAP) experiment [4] suggests that about 85% of the 
total matter content of the universe is dark. This constitutes 23% of the total content of the 
universe. The rest 73% is the dark energy, whereas the remaining 4% is the known luminous 
matter. Nature and identity of the constituents of this non-luminous matter is mostly unknown. 
However, the indirect evidences suggest that most of them are stable , nonrelativistic (Cold Dark 
Matter or CDM) and Weakly Interacting Massive Particles (WIMPs) [5, 6, 7, 8]. Despite the 
wide success spectrum of Standard Model (SM), explanation of CDM poses a challenge to SM 
as no viable candidate for CDMs has been obtained within the framework of SM. Hence, there 
are attempts to explain the DM WIMP candidates with theories beyond standard model. 

Phenomenology of different extensions of scalar sector of the SM had been explored by many 
groups [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. However, addition of one real scalar singlet 
to the SM provides the simplest possible minimal renormalizable extension to the scalar sector 
of SM. In addition, invoking a Zi discrete symmetry under which the additional singlet is odd, 
gives rise to the singlet as a viable DM candidate. In this work, we explore the parameter space 
of the model to accommodate the results of different dark matter direct detection experiments. 
We further predict the event rates and the annual variation of the event rates to be observed 
by the Liquid Argon Detector experiment. In the experiments for direct detection of DM, the 
WIMP scatters off the target nucleus of the material of the detector giving rise to recoil of the 
nucleus. The energy of this nuclear recoil is very low (~ keV). The signal generated by the 
nuclear recoil is measured for direct detection of dark matter. 

There are several ongoing experiments for direct dark matter searches. Some of them are 
cryogenic detectors where the detector material such as Germanium are kept at a very low 
temperature background and the nuclear recoil energy is measured using scintillation, phonon 
or ionization techniques. The experiments like CDMS (Cryogenic Dark Matter Search uses 
Germanium as detector material) at Soudan Mine, Minnesota [27] use both ionization and phonon 
techniques. In phonon technique, the energy of the recoil nucleus sets up a vibration of the 
detector material (Ge crystal for CDMS). These vibrations or phonons propagate at the surface 
of the detector crystal and excites quasi-particle states at materials used in the pulse pick up 
device. Finally the heat produced by these quasi-particle states is converted to pulses by SQUID 
(Superconducting Quantum Interference Device) amplifiers. CDMS carries out two experiments 
- one with Germanium and the other with Silicon in order to separate the neutron background. 
As Germanium nucleus is heavier (A = 73) than Silicon (A = 28), WIMPs interact with 73 Ge 
with higher probability than with 28 Si, but neutron being strongly interacting will not make any 
such discrimination. Thus any excess signal at the 73 Ge detector over the 28 Si in CDMS will be 
a possible signature for dark matter. 

The DAMA experiment at Gran Sasso (uses diatomic Nal as the detector material) [28], 
uses the scintillation technique for detection of the recoil energy. There are other class of liquid 
or gas (generally noble gases) detectors that measure the recoil energy by the ionization of the 
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detector gas. The ionization yield is amplified by an avalanche process and the drifting of these 
charge reaches the top (along z-axis) where they are collected by the electrodes for generating 
a signal. These types of detectors known as TPC (Time Projection Chamber) are gaining lot 
of interest in present time for their better effectiveness and resolution in detecting such direct 
signals of recoil energy from a DM-nucleon scattering. As mentioned, they generally use noble 
gases such as Xenon, Argon or Fluorine etc. The XENON- 10 experiment [29] at Gran Sasso is a 
liquid Xenon TPC with target mass of 13.7 Kg, whereas its upgraded version, the XENON-100 
experiment [30] uses 100 Kg of the target mass. The CoGeNT (Coherent Germanium Neutrino 
Technology) experiment [31] also uses Ge as detector material and designed to detect dark matter 
particles with masses less than that to be probed in CDMS. Recently the CoGeNT collaboration 
has reported an excess of events above the expected background [31]. The experiment ArDM 
(Argon Dark Matter Experiment) [32] plans to use 1 ton of liquid 39 Ar gas for the TPC. There 
are other experiments that use other techniques like PICASSO [33] etc. at SNOlab in Canada 
but here we consider CDMS, DAMA, CoGeNT and Xenon experiments for the present study. 
We restrict the relevant couplings of the scalar dark matter by using the bounds on dark matter- 
nucleon scattering cross sections from these three experiments and further, utilising the WMAP 
experimental data. We predict possible direct detection event rate as well as annual variation of 
event rate in liquid ArDM experiment for different possible dark matter masses. 

The paper is organized as follows. In Section 2 we briefly discuss the model. The CDMS, 
XENON, CoGeNT and DAMA bounds used to constrain the parameter space of the model is 
described in Section 3. In Section 4 we present the formalism of Direct detection rate calculations 
and compute such rates for liquid Ar detector. Section 5 contains summary and conclusions. 



2 Singlet extended Standard Model : A brief outline 

The framework of the simplest scalar sector extension of the SM involving addition of a real 
scalar singlet field to the SM Lagrangian has been discussed in detail in [24, 25]. In this section 
we present a brief outline of the model and emphasize on those of its aspects that would be 
relevant for discussions to follow. 

The most general form of the potential appearing in the Lagrangian density for scalar sector 
of this model is given by 1 

V(H,S) = —H^H+-(H^H) 2 + 6 -±H^HS+-H^HS 2 
2 4 2 2 

+ ^ + ^ (i) 

where H is the complex Higgs field (an SU (2) doublet) and S is a real scalar gauge singlet that 
defines our minimal extension to the scalar sector of SM. The singlet S needs to be stable in 
order to be considered as a viable dark matter candidate. Stability of S is achieved within the 
theoretical framework of the model by assuming the potential to exhibit a Z 2 symmetry S — > —S. 




1 We used same notations as used in [24, 25]. 
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Figure 1: Diagram for singlet-nucleon elastic scattering via higgs mediation 



This ensures absence of vertices involving odd number of singlet fields S (Si = K3 =0). Using 
unitary gauge, we define 



H = (^J (2) 

where h is the physical Higgs field and v = 246 GeV is the VEV of the H scalar determined by 
the parameters m and A as v = \J ~ 2 ™ 2 ■ The mass terms of the two scalar fields h and S are 
identified as 

V mass = l -(M 2 h h 2 + M 2 S S 2 ) (3) 

where, 

M 2 h = -m 2 = \v 2 /2 

Ml = k 2 + S 2 v 2 /2 (A) 

The scalar field S is stable as long as the Z 2 symmetry is unbroken and appears to be a candidate 
for cold dark matter in the universe. 

In the present work, we investigate the prospect of such a candidate in direct detection 
experiments through its scattering off nucleon (N) relevant for the detector. The lowest order 
diagram for the process has been shown in Fig. 1. The cross section corresponding to the elastic 
scattering (SN — > SN) in the non-relativistic limit is given by [11] 

^ = ^Mn\ 2 ( ml \ 
N 4tt \M S 2 M h V 1 ^ 

where, m r (N, S) = MnMs/(Mn + Ms) is the reduced mass for the target nucleus of the two body 
scattering .SW — > SN and A n is the relevant matrix element. For the case of non-relativistic 
nucleons the singlet-nucleus and singlet-nucleon elastic scattering cross sections are related by 

[11] 

scalar A 2 m 2 r (nudeUS, S) scalar 

nucleus - w? r (nuc\eon, S) nuclcon K ) 
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Figure 2: Dashed lines: Scalar singlet-nucleon elastic scattering cross section (f"nudeon) ^ ^ 
function of the singlet mass for different values of [ <5s | with 2 different values of higgs mass Mh, 
114 GeV(left panel) and 200 GeV(right panel). Solid lines: 90% CL. experimental upper limits 
on a™ leon from CDMS 2009 Ge, CDMS Soudan(all), XENON-10 and XENON-100. The shaded 
areas corresponding to 99% C.L. regions allowed from DAMA (with and without channelling) 
has also been shown. 



Numerically evaluating the matrix element appearing in Eq. (5) the singlet-nucleon elastic scat- 
tering cross section can be written as [11] 

scalar , R , 2 ( 100_GeV\ V50GeV\ 2 1n -42 2^ m 

It is evident from the above equations that the scalar cross section for scalar singlet depends on 
two couplings 5 2 and k 2 . Hence the dark matter direct detection rate with scalar singlet as the 
dark matter candidate also depends on the two couplings 5 2 and k 2 . The constraint on direct 
detection rates of dark matter from different experiments thus can also put constraint on the 
parameter space (5 2 , k 2 ) of scalar singlet dark matter. In the next section we will discuss how 
some of the recent direct detection experiments of dark matter can put limit on the scalar singlet 
model of dark matter. 



3 Constraining the model parameters 

The direct detection rates of scalar singlet dark matter is governed by the scalar-nucleon cross 
section (a^ieon) f° r a given scalar singlet mass M s . Direct detection experiments like CDMS, 
DAMA, Xenon thus can constraint the scalar singlet parameter space. Preliminary studies of 
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scalar singlet dark matter concerning CDMS, Xenon are given in Refs. [12, 18, 19]. Also, 
explanation of DAMA bounds with scalar singlet dark matter is described in [34, 20]. In this 
work we make a detailed study of all the direct detection bounds together with the relic density 
limits from WMAP and find the relevant parameter space for the scalar singlet to be a successful 
dark matter candidate. 

3.1 Constraints on 62 as a function of dark matter mass 

The coupling 82 is of severe importance as this is the only coupling in this model which determines 
the annihilation of the scalar to other standard particles. Moreover given the singlet mass (Ms), 
8 2 is the sole coupling that controls the singlet-nucleon cross section and the quadratic dependence 
(Eq. 7 ) in particular reflects that this cross section is insensitive to the sign of 8 2 . 

In Fig. 2 the singlet-nucleon elastic scattering cross section is plotted as a function of singlet 
mass (dark matter mass) Ms for different values of 82- The plots are presented at two different 
values of Higgs mass (M^), namely 114 GeV (left panel) and 200 GeV (right panel) of Fig. 2. For 
comparison, the 90% C.L. (confidence limit) results obtained from CDMS II experiment (CDMS 
2009 (Ge)) [26] are plotted in Fig. 2 with the similar results from the combined analysis of full 
data set of Soudan CDMS II results (CDMS Soudan (All)) [26], XENON-10 [29] and XENON- 
100 experiment [30]. One sees from Fig. 2 that the DAMA and CoGeNT results are constrained 
in closed allowed regions unlike the other experiments that provide upper bounds of the allowed 
masses and cross-sections of dark matter. In this context it may be noted that the consideration 
of ion channelling in Nal crystal in DAMA experiment is crucial in the interpretation of its results. 
The presence of channelling in Nal, affects the allowed mass and cross section regions of the DM 
particles inferred from the observation of an annual modulation by the DAMA collaboration. 
The effect of channelling has been discussed extensively in [35]. For our analysis we consider the 
allowed mass-cross section limits inferred from observed annual modulation of DAMA for both 
the cases — with channelling and without channelling. From Fig. 2 it is seen that for larger values 
of 8 2 , higher singlet mass domain is required in order to represent the experimentally allowed 
region for ^dcon - DM mass plane. Also, the region of overlap of cr^teon — M s plots for a fixed 
value of 82, becomes larger for higher Higgs mass. 

In Fig. 3 we present plots for upper limits of the range of \8 2 \ (as a function of singlet mass) 
that would reproduce cross section values (computed with Eq. 7) below the 90% CL limits of 
different experiments shown in Fig. 2. It is evident from the plots that, as in Fig. 2, the allowed 
values of the coupling 82 is sensitive towards the Higgs mass. Moreover the minimum allowed 
value of dark matter mass in this model is also dependent on Higgs mass. The appearance of 
local minima at low Ms domain of the plots are due to the behaviour of experimental bounds 
that show a sudden upturn followed by a rise of the curves at low M s (M s — > 0) region. For 
Higgs mass of 120 GeV (114 GeV) this minima corresponds to the \8 2 \ value ~ 0.5(0.1). 

3.2 Constraining the 62 — K2 parameter space 

In the previous subsection, constraining of the parameter 8 2 using other direct detection experi- 
mental results are addressed. But as described earlier, the present dark matter model with scalar 



6 



2.5 - 



I 

CDMS2009(Ge) 
CDMS Soudan (All) 
XENON- 10 
XENON- 100 



iO 1.5- 



0.5 - 



200 300 

M s (in GeV) 



M h = 114GeV 

n 



2.5 



2t 



to 1.5 



0.5 



■■■ CDMS 2009 (Ge) / 
■- CDMS Soudan (All) / 
— XENON-10 
-- XENON-100 



• ; 
i 

I r 

I I 



i ■ I / 
' l 
i ■ i I 

!/// 



f A 
.•* * 
Is' 



100 



200 



(in GeV) 



300 



M, = 200 GeV 



411(1 



501) 



Figure 3: Upper limits on | <5^2 1 as a function of dark matter mass from 90% CL experimental 
bounds on spin-independent WIMP-nucleon cross section. Plots are done for 2 values of Higgs 
mass: 114 GeV (left panel) and 200 GeV (right panel) . 



singlet also depends on the other parameter namely k 2 (Eq. (4)). From Eq. (4) one sees that 
the singlet mass term Ms is governed by two model parameters namely <5 2 and k 2 . For a given 
value of 62, real values of Ms can be obtained only by excluding all k 2 values less than — <5 2 -y 2 /2. 
In principle any value of dark matter mass Ms can accommodate all (62, K2) values satisfying 
K2 + 82V 2 / 2 = Ms 2 . However at large K2 values with K2 ^> c^f 2 , the singlet mass is predomi- 
nantly K2 driven and is scaled with it as ^/K2■ The interplay between 62 and K2 in setting a given 
singlet mass is represented in Fig. 4 where we plotted (dashed lines) different iso-M^ contours in 
82 — sign(/t 2 )|/t2| 1 ^ 2 plane. For a given singlet mass Mg, the range of | <^2 1 consistent with CDMS 
II/XENON-10 limits on WIMP-nucleon scattering cross section (as discussed in Sec. 3.1) would 
then correspond to a segment of the corresponding iso-M^ contour in 82 — sign^)!^! 1 / 2 plane. 
Its projection on the k 2 axis would give the corresponding range of the parameter k 2 . In Fig. 
4 the shaded regions represent the domains of the model-parameter space <5 2 — sign(K 2 )|K 2 | 1//2 , 
that is consistent with 90% CL limits on WIMP-nucleon elastic scattering cross sections from 
analysis of CDMS Soudan (All) (left panel), XENON-10 (middle panel) and XENON-100 (right 
panel) results. The dependence of the allowed model-parameter space on Higgs mass is shown 
by plotting the allowed areas for three different values of Higgs mass namely 114 GeV, 150 GeV 
and 200 GeV. The region of the parameter space to the left of Ms = line is excluded as points 
(82, k 2 ) in that region would give negative values of mass. Moving from lower Ms to higher Ms 
domain in the parameter space allows more and more room for (<5 2 , k 2 ) to represent DM-nucleon 
scattering cross section consistent with its experimental bounds - a feature also apparent from 
Fig. 2 and Fig. 3 (discussed in Sec. 3.1). At very low scalar singlet mass regime (Ms ~ — 10 
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Figure 4: Shaded region: Range of the parameter space 8% — sign(fc 2 )|ft2| 1/ ' 2 consistent with 
90% CL limits of WIMP-nucleon scattering cross section from CDMS Soudan (All) (left panel), 
XENON- 10 (middle panel) and XENON- 100 (right panel) corresponding to three different values 
of Higgs mass - 114 GeV, 150GeV, and 200 GeV. Dashed lines: Different iso-Ms contours in 
§2— sign^)!^! 1 / 2 plane. The region spanned by red dots describe the values of model parameters 
consistent with WMAP measurements of dark matter relic density : 0.099 < Qh 2 < 0.123. 



GeV) higher peaked value of the WIMP-nucleon scattering cross section-limit concedes a thread 
like extension of the allowed parameter space along the corresponding iso-M5 contours. The 
sudden drop of the experimental limits of WIMP-nucleon cross sections in 10 — 20 GeV mass 
regime severely restricts the width of the above mentioned thread-like extension of the allowed 
parameter space at Ms ~ — 10 GeV. This drop is more robust for the case of the experiment 
denoted in this work as "CDMS Soudan (All)" than that of "XENON" experiment leading to a 
more prominent appearance of the thread-like zone in left panel of Fig. 4. The 82 — ^2 parameter 
zone obtained from WMAP constraint [24] (0.099 < Qh 2 < 0.123, Q being the dark matter relic 
density and h is the Hubble parameter normalized to 100 Km sec -1 Mpc -1 ) are also shown by 
the red colored dots in Fig. 4, for comparison. From Fig. 4 it is seen that CDMS/XENON 
upper limits of WIMP-nucleon scattering cross section together with WMAP observation allows 
only small \8z\ regime (^ 0.2) of the model parameter space for different higher values of Ms , 
although for very small Ms values (0 — 10 GeV) CDMS limit concedes more room for {82 | (upto 
~ 1.0) 

The results of DAMA experiment restricts the variations of cr^icon with DM mass in two small 
contours represented as shaded areas in Fig. 2. They represent the 99% C.L. regions in (c^^, n - 
DM mass) space. Interpretation of DAMA results with channelling (without channelling) requires 
WIMP-nucleon scattering cross sections of order 10 _41 cm 2 (10~ 40 cm 2 ) along with two locally 
preferred zones of dark matter mass - one around ~ 12 GeV (referred as DAMA-a in this work) 
and the other around ~ 70 GeV (referred as DAMA-b in this work). The DAMA solution 
corresponding to the large DM mass regime (DAMA-b) is completely excluded by observed 
limits on the WIMP-nucleon scattering cross section form other direct detection experiments 
like CDMS, XENON etc. The other DAMA solutions corresponding to lower DM mass regime 
(DAMA-a) are also largely disfavoured by XENON. However though a small region at the lower 
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Figure 5: Left panel: (Shaded region) Range of the parameter space 82 — sign(K 2 )|K 2 | 1//2 consistent 
with 90% CL limits of WIMP-nucleon scattering cross section from CDMS 2009 (Ge) for three 
different values of Higgs mass 114 GeV, 150 GeV and 200 GeV. Different iso-Ms contours in 
82 — sign(/=c 2 ) 1^2 1 1,/2 plane are shown by dashed lines. The region spanned by red dots describe the 
values of model parameters consistent with WMAP measurements of dark matter relic density : 
0.099 < Qh 2 < 0.123. Middle panel: (small thread-like regions outlined by dashed lines) Region 
in parameter space consistent with (CDMS 2009 Ge + DAMA(with channelling)) limits for 
= 114 GeV (1 st column), 150 GeV (2 nd column), 200 GeV (3 rd column). The region allowed 
from only CDMS 2009 (Ge) limits are shaded for corresponding values of Higgs mass. The 
contour for Ms = 10 GeV are shown by dotted lines. Right panel: Range of the parameter space 
62 — sign(ft2)|ft2| 1//2 consistent with 90% CL limits of WIMP-nucleon scattering corresponding 
to the lower DM mass regime allowed from DAMA (without channelling) (upper panel) and 
CoGeNT (lower panel) for three different values of Higgs mass 114 GeV, 150 GeV and 200 GeV 



DM mass regime ("with channelling" case) in DAMA-a zone is barely consistent with 90% CL. 
with CDMS experiment denoted in this work as "CDMS 2009 (Ge)", the corresponding limit 
"CDMS Soudan (All)" experiment excludes the entire DAMA-a zone. In the left panel of Fig. 5 
we have shown region of 82 — sign^)]^! 1 ^ 2 parameter space that corresponds to singlet-nucleon 
elastic scattering cross sections within its 90% CL. limit from CDMS 2009 (Ge). The middle 
panel of Fig. 5 shows the 82 — sign^)!^! 1 / 2 parameter space, consistent with both CDMS 2009 
(Ge) and DAMA results (with channelling). The small thread like regions (marked red and 
indicated within closed dashed lines) in the middle panel of Fig. 5 represent the parameter space 
domain that fits singlet-nucleon scattering cross section with DM-nucleon cross section consistent 
with both CDMS 2009 (Ge) and DAMA limits (with channelling). These are presented in three 
different columns that correspond to three values of Higgs mass namely 114 GeV, 150 GeV and 
200 GeV. The parameter space consistent with CDMS 2009 (Ge) limit only for corresponding 
values of Higgs masses are also shown in the respective columns by grey shades. The iso-Ms 
contour for M$ = 10 GeV, spanning through the DAMA + CDMS 2009 (Ge) allowed regimes, is 
also shown in each column of the same figure to illustrate the fact that parameter space regions 
consistent with both DAMA and CDMS 2009 (Ge) correspond to a singlet dark matter mass 
of around ~ 10 GeV. In the right panel of Fig. 5 we have shown the allowed region in the 
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52 — sign(A€ 2 ) | 1 1/72 parameter space consistent with DAMA-a region without channelling (upper 
panel) and the CoGeNT result (the lower panel). 

4 Predictions for Rates at Argon detector 

The Argon detector is a Noble liquid detector where the liquefied noble gas Argon is used as 
targets for direct detection of WIMPs. Because of the high density and high atomic number the 
event rate is expected to be large. Also because of its high scintillation and ionization yields 
owing to its low ionization potentials, it can effectively discriminate nuclear recoils and other 
backgrounds from 7 or electrons. The ArDM (Argon Dark Matter) experiment at surface of 
CERN uses one such detector that envisages one ton of liquid Argon in a cylindrical container. 
This has a provision for three dimensional imaging for every event. A strong electric field along 
the axis of the cylinder helps drifting of the charge - produced due to the ionization of liquid 
Argon by WIMP induced nuclear recoil - to the surface of the liquid. This charge then enters 
into the gaseous phase of the detector (TPC) where it is multiplied through avalanche and finally 
recorded by a position sensitive readout. 

In this section we estimate WIMP signal rates for such Argon detector. The differential rate 
for dark matter scattering detected per unit detector mass can be written as 

dR , r ^ da f „, . , 

n t®^t^ / f(v)dv (8) 



d\q\ 2 d\q\' 

where is number of target nuclei per unit mass of the detector, $ is the dark matter flux, f(v) 
denotes the distribution of dark matter velocity v (in earth's frame). The integration is over all 
possible kinematic configurations in the scattering process. |q| is the momentum transferred to 
the nucleus in dark matter-nucleus scattering and a being the corresponding cross section. The 
recoil energy of the scattered nucleus can be expressed in terms of the momentum transfer |q| as 

E R = |q| 2 /2Mjv = m 2 r v 2 (l -cos9)/M N (9) 

where M N is the nuclear mass, 9 is the scattering angle in dark matter - nucleus center of 
momentum frame and m r is the reduced mass given by 

M N M S 

m r = T7 7T- ■ (1°) 

M N + M S v 1 

In the above M$ is the dark matter mass. Expressing the dark matter flux $ in terms of the 
local dark matter density p s , velocity v and mass M s . With N T = l/M N and writing |q| 2 in 
terms of nuclear recoil energy E R , Eq. (8) can be rewritten as 

2 ,, , 2 / vf(v)dv (11) 



dE R M s d\q\ 
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where 



~M N E R 
2m 2 . 



1/2 



(12) 



The dark matter - nucleus differential cross-section for the scalar interaction is given by [5] 



da 



a 



scalar 



<Z|q 



2 4m 2 ed v 2 



F 2 {E R ) . 



(13) 



Here cr scalar is dark matter-nucleus scalar cross-section and F(E R ) is nuclear form factor given 
by [36, 37] 



F(E R 



3ji(qRi) 



cxp 



' q 2 s 2 ^ 



(14) 



qRi 

Rl = {r 2 -hs 2 fl 2 
r = 1.2A 1 ' 3 

where s(~ 1 fm) is the thickness parameter of the nuclear surface, A is the mass number of the 
nucleus, ji(qRi) is the spherical Bessel function of index 1 and q — |q| = \f2M^E R as from Eq. 
(9). Assuming distribution /(v ga i) of dark matter velocities (%a) with respect to galactic rest 
frame to be Maxwellian, one can obtain the distribution f(v) of dark matter velocity (v) with 
respect to earth rest frame by making the transformation 



v = v 



gal ~~ v e 



(15) 



where v e is the velocity of earth with respect to Galactic rest frame and is given as a function 
of time t by 



Vq + V orb COS 7 COS 



^ 27T(t-t ) ' 



(16) 



In Eq. (16) v Q is the speed of the solar system in galactic rest frame, T (1 year) us the period of 
earth's rotation about sun , t = 2 nd June (the time of the year when the orbital velocity of earth 
and velocity of solar system point in the same direction) and 7 ~ 60° is the angle subtended by 
earth orbital plane at Galactic plane. The speed of solar system v Q in the Galactic rest frame is 
given by 



Vq = V + V- 



pec 



(17) 



Vq being the circular velocity of the local system at the position of solar system and v pec = 
12km/sec, called peculiar velocity, is speed of solar system with respect to the local system. 
Physical range of v is 170 km/sec < v < 270 km/sec (90 % C.L.) [38, 39]. In this work we 
consider the central value - 220 km/sec for v . The term cos2[-7r(t — t )/T] in the velocity is 
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responsible for annual modulation of dark matter signal. Introducing a dimensionless quantity 
T(E R ) as 



T(E R ) = 



which can also be expressed as 



T{E R ) 



Ave 



-Vq 



erf 



A jv 



+ Vn 



f(v) 



we obtain from Eqs. (11) and (13) 



dR 



dE 



scalar 

PS :FHE R ) 



R 



A V(B Msm 2 r ' 



erf 



— erf 



dv 



(18) 



+ 



Vq 



erf 



(19) 



(20) 



The local dark matter density p s may be taken as 0.3 GeV/cm 3 . The observed recoil energy 
(E) in the measured response of the detector is a fraction (Qx) of actual recoil energy (E R ) at 
the time of scattering. This fraction Qx = EjE R (called as quenching factor) is different for 
different scattered nucleus X. For 39 Ar, Q^ r = 0.76. Thus the differential rate in terms of the 
observed recoil energy E for 39 Ar detector can be expressed as 



AR 
~AE 



(E) 



(E + AE)/Q Al d R A AE R 

e/q Ai dE R [ R) AE 



(21) 



In left panel of Fig. 6 we show the expected differential rates (/kg/ day/keV) for different observed 
recoil energies in Argon detector considering scalar singlet as the dark matter candidate. Four 
representative cases have been plotted and they are denoted as (a), (b), (c) and (d). For all the 
plots, (a) - (d), the chosen values of coupling 62 (as also corresponding scalar singlet masses, Ms) 
are consistent with current CDMS and WMAP limits. All the plots show that the rate falls off 
with the increase of recoil energy. Plots (a) and (b) are the variations of rates for the same set 
of Higgs mass (Mjj = 114 GeV) and singlet mass (Ms = 100 GeV) but for different values of the 
coupling 5 2 (S 2 = — 0.1(— 0.03) for plot (a) ((b)). Plots (a) and (b) show a decrease of the rate 
when \5 2 \ decreases. For example, in case of recoil energy E = 50 GeV the calculated rates from 
plots (a) and (b) are 8.5 x 10~ 6 (for \5 2 \ = 0.1) and 7.2 x 10~ 7 (for \5 2 \ = 0.03 ) respectively in the 
units of /kg/day/keV. Plots (c) and (d) compare the variation of rates for two different scalar 
masses namely Ms = 200 GeV (plot (c)) and Ms = 100 GeV (plot (d)) for same values of M^, 5 2 
(200 GeV, —0.03). From plots (c) and (d) it is seen that the rate increases for any particular 
value of recoil energy with decrease of dark matter mass. For example in case of E = 50 GeV, 
the rates for M s = 100 GeV (plot c) and M s = 200 GeV (plot d) the calculated rates are 8 x 10~ 8 
and 1.4 x 10~ 8 respectively in the units of /kg/day/keV. This is evident from the expression for 
scalar cross section (Eq. (7) which varies as M$ 2 and direct detection rates is linear with the 
scalar cross section. One can compare plots (b) and (d) to see the effect of Higgs mass values 
on the rate. For Ms, S 2 (100 GeV, —0.03) , the estimated rates in the present calculations at 
E = 50 keV are 7.2 x 10~ 7 (/kg/day/keV) for M h = 114 GeV (plot (b)) which is reduced to 
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E (keV) Days in a year 



Figure 6: Left Panel: Plot of predictions for dark matter detection rates (per kg per day per keV) 
in Argon detector as a function of observed recoil energy. The plots are shown for two values 
of Mh - 114 GeV and 200 GeV and for two sets of (Ms, 62) values which are consistent with 
CDMS limits as well as observed relic density of dark matter (WMAP). In the inset we show 
the corresponding plot for M h = 114 GeV and for two sets of (Ms, 62) values simultaneously 
consistent with limits on scattering cross section from CDMS 2009 (Ge) and DAMA and also 
with WMAP. Right panel: Predicted annual variation of event rates in Argon detector over one 
year for Mh = 114 (GeV). The upper panel corresponds to a (Ms, 62) value consistent with 
CDMS+DAMA+WMAP while the lower panel corresponds to a (Ms, 5 2 ) value consistent with 
CDMS+WMAP. 

1.4 x 10 -8 (/kg/day/keV) for Mh = 200 GeV (plot (d)). In the inset of this same figure we show 
the calculated prediction of rates for Ms = 10 GeV, 62 = 0.4. This is compatible with the CDMS 
and DAMA bound together, other than satisfying the WMAP limits. The Higgs mass is kept at 
114 GeV. As an example, for E = 30 GeV the calculated expected event is 1.4 x 10 -3 per day 
for a ton of the detector. 

One very positive signature of dark matter in direct detection method is the annual variation 
of detection rate. This periodic variation arises due to periodic motion of the earth about the sun 
in which the directionality of earth's motion changes continually over the year. As a result, there 
is an annual variation in the amount of dark matter encountered by the earth. The detection 
of annual variation in direct detection experiments serve as a smoking gun signal for existence 
of dark matter. In right panel of Fig. 6 we show the calculated annual variation of event rate 
(/kg/day) in different times of a year. In the upper panel, we have chosen the parameter set 
that is compatible with CDMS, DAMA and WMAP limits (as in inset of the left panel) and for 
the lower panel we have given the results for M h = 114 GeV, M s = 100 GeV and 5 2 = —0.03. 
As expected, the yield is maximum on 2nd June when the direction of motion of the earth is the 
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same as that of the solar system. 



5 Summary and Conclusions 

In the present work we consider a simplest extension of SM, introducing a real scalar singlet along 
with a discrete Z 2 symmetry which ensures stability of the singlet. Such singlet is considered 
as a viable cold dark matter candidate. The scattering of this singlet dark matter off nuclei 
in the detector can be observed by measuring the energy of the recoil nuclei. The calculated 
singlet-nucleon scattering cross section in this model explicitly depends on the coupling 5 2 and 
implicitly on k 2 (as defined in Eq. (1)). We constrain the 5 2 — ^2 parameter space using the 
recent bounds on the WIMP-nucleon scalar cross section as function of WIMP mass, reported by 
the CDMS collaboration and also reported by the XENON-10 and XENON-100 collaborations. 
The allowed zones in the parameter space follow a typical pattern determined by the shape of 
the WIMP mass dependence of the experimental limits considered here and the way the scalar 
singlet mass is related to 5 2 and k 2 . The allowed zones vary for different Higgs masses and they 
are consistent with WMAP limits. 

We investigate the effect of the inclusion of DAMA results (with channelling) on the 5 2 — k 2 
parameter space but the allowed zone is found to be extremely small and representative of a 
very low dark matter mass (around ~ 10 GeV). We also compute the range of 5 2 — k 2 parameter 
space consistent with DAMA (without channelling) and CoGeNT experiments. 

Utilising the constrained parameter space we estimate the possible detection rates and their 
annual variations for a liquid Argon detector. 
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